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Abstract 

Traditional derivation of Gibbs canonical distribution and the justification 
of thermodynamics are based on the assumption concerning an isoenergetic 
ergodicity of a system of n weakly interacting identical subsystems and pas- 
sage to the limit n — > oo. In the presented work we develop another approach 
to these problems assuming that n is fixed and n>2. The ergodic hypothesis 
(which frequently is not valid due to known results of the KAM-theory) is 
substituted by a weaker assumption that the perturbed system does not have 
additional first integrals independent of the energy integral. The proof of 
nonintegrability of perturbed Hamiltonian systems is based on the Poincare 
method. Moreover, we use the natural Gibbs assumption concerning a ther- 
modynamic equilibrium of subsystems at vanishing interaction. The general 
results are applied to the system of the weakly connected pendula. The aver- 
aging with respect to the Gibbs measure allows to pass from usual dynamics 
of mechanical systems to the classical thermodynamic model. 



1 Introduction 

The classical approach to the justification of thermodynamics is based on the use of 
Gibbs canonical distribution 

e -0H 

P ^ V)= fe-Wdxdy - (L1) 

Here x and y are canonical coordinates and momenta respectively, H(x, y) is the 
Hamiltonian function of a mechanical system, (3 = const. The function p is treated 
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as a stationary density of probability distribution. More precisely, we suppose that 
the system in the given state is a random event, and the probability of detection of 
the system in area D of the phase space is equal to 

J p{x,y)dx dy. 

D 

The parameter j3 is determined from the equality 

E = J Hpdxdy, (1.2) 

where E is average energy of the system. Usually this parameter is supposed to be 
(3 = Ty, where T is the absolute temperature, k is the Boltzmann constant PQ [I]. 
We assume that the integrals (1.1) and (1.2), extended to the whole phase space, 
are converging. 

The justification of Gibbs distribution is one of basic problems of statistical 
mechanics. For this purpose one usually considers N indistinguishable systems with 
the Hamiltonians 

H(xM,yM), a = l,...,N, (1.3) 

where x^ a \ y^ (1 < % < n) are canonical coordinates of a system with number a. 
Let us emphasize that n and H do not depend on a. 

Let us introduce the unified phase space of dimension 2nN with canonical vari- 

^ (X ) y) = (x«yW..,xW^) 

and the Hamiltonian 

JV 

H £ (X,Y) = Y,H{x^\y^)+eH x {X,Y,e) , (1.4) 

a=l 

where e is a small parameter, which then will approach to zero. The system with 
the Hamiltonian (1.4) describes the dynamics of N weakly interacting subsystems 
with the Hamiltonians (1.3), and the function eH\ denotes the interaction energy 
of subsystems. 

The basic idea of derivation of Gibbs distribution is the assumption concerning 
the ergodicity of the system with the Hamiltonian (1.4) at e ^ on a level surface 
H £ = NE [3]. Let / be an integrable function on a phase space of a system with 
n degrees of freedom. According to an indistinguishability principle, in statistical 
mechanics one considers only such functions of X and Y, which do not vary under 
permutations of groups of variables x' a \ y^ . Let us assume that 

nx,Y) = ^Y,^ x(a) ^ {a) )- ( L5 ) 



2 



Such functions are frequently called summators. Using ergodic hypothesis, we obtain 
for theaverage with respect to time (1.5) in the limit as e — > the formula 

(F) = J f(x,y)pN(x,y)dxdy. 

The explicit expression for the density pN can be found in t 3j. With some additional 
suppositions Z\ we can prove that p n — > p as N — > oo and for the average of function 
/ along solutions of the system with the Hamiltonian (1.3) the formula 

(/) = j fpdxdy 

is valid. 

On this way one has a series of fundamental difficulties. The main difficulty 
is the justification of the ergodic hypothesis. Only recently the ergodicity of some 
simplified models has been established jS]. Since the Gibbs distribution (1.1) does 
not depend on the form of an interaction energy H\, F.A.Berezin [3] stated the 
idea of weakening of the ergodic hypothesis: it is enough to require ergodicity of 
the system with the Hamiltonian (1.4) for e ^ for set of functions H\ everywhere 
dense. However, if the interaction energy H% has no singularities and surfaces of the 
energy level 7i £ = NE are compact, the ergodic hypothesis (even in the weakened 
variant) isrefuted by the KAM-theory 

EXAMPLE 1. Let us consider a system with the Hamiltonian, following [3| 



n 



\ Y,{£ + "Ye*) + eV*{x) , (1.6) 



where u = const ^ 0, V4 > is a polynomial of the 4-th degree of coordinates 
xi, . . . , xjsr, e > 0. Such systems play an essential role in the theory of heat capacities 
of rigid bodies [2]. The coefficients of non-negative polynomials of the 4-th degree 
form some set A in a finite-dimensional space of coefficients of all polynomials of 
the 4-th degree from iV variables. F. A. Berezin [S] has raised the following question: 
Is it correct, that for almost all points A the system with the Hamiltonian (1.6) is 
ergodic? The answer is negative. Indeed, let us consider polynomials of the form 



V4 = a a Xa + W&(x) , a a > 



and let the coefficients of the polynomial W4 be small. If W4 = 0, the system with 
the Hamiltonian (1.6) is integrated by a separation of variables. If we change to 
action-angle variables /„, ip a in a system with one degree of freedom 

_ (x 2 + uj 2 y 2 ) 4 



then dh d 2 h 

*» = Sr>o> x * = ^w>°- ( L7 ) 
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In the last inequality we essentially use the supposition, that ea a > 0. 
Let H = J2K{Ia)- If 



d 2 H 



dl 



# o 



^0 



'U 



then, according to the KAM-theory jH], the system with the Hamiltonian (1.6) is 
not ergodic on each positive energy level for small W± (when coefficients of W4 are 
located in some small neighborhood of zero): the most part of this surface is foliated 
on iV-dimensional invariant tori with conditional-periodic trajectories. In our case 
the determinant (1.8) is equal to 



— Ai . . . A n 



Ai 



+ . . . + 



An 



what does not equal zero on account of (1.7). 

It is still necessary to add that according to the above mentioned approach to 
the derivation of the Gibbs distribution (due to Darwin and Fouler (see [3])) the 
expression for the density of probability distribution (1.1) does not depend on the 
interaction of subsystems. While such an interaction is always present. 

In this work we develop another approach to the justification of the formula for 
probability density due to classical works by Gibbs pp. This approach will be based 
on the theory of integrability of Hamiltonian systems [7] . 



2 The main result 

The statistical approach to the theory of dynamical systems assumes a refusal of 
consideration of separate trajectories. Instead, in the phase space of a system 

x = v{x) (2.1) 

we introduce the density of probability distribution p(x, t) > 0, explicitly depending 
on time in general case. Let D be any measurable area of the phase space, g l be a 
phase flow of the system (2.1). This flow can be presented as a stationary flow of a 
fluid. Since the area g t {D) consists of the same moving phase points, it is natural 
to assume that the probability of detection of the system in the area g t {D) does not 
depend on t. Hence, this probability 

J p(x,t)dx 

9*(D) 

will be an integral invariant of the system (2.1). But then the density of probability 
distribution satisfies the Liouville equation: 

dp 

— + div(pv) = . 
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If divf = 0, then p is a first integral of the system (2.1). This condition, 
obviously, holds for the Hamiltonian systems of differential equations. It is natural 
to consider stationary distributions for autonomous systems, when p does not depend 
on t explicitly. The detailed discussion of these problems can be found in [21 II]- 

Locally, in a small neighbourhood of a nonsingular point the dynamical system 
(2.1) possesses a complete set of independent first integrals (in amount of m — 1, 
where m is a dimension of the phase space). However, in a typical situation they 
can not be extended to single- valued integrals defined in the whole phase space. The 
density of probability distribution is a single-valued function on the phase space by 
definition. Generally Hamiltonian systems with compact level surfaces do not admit 
integrals independent of the Hamiltonian H, though they are not ergodic (on energy 
levels) [7j. This argument allows to state at once in some important cases, that 

P = /(#), (2-2) 

and to reduce the problem to determination of the form of function /. 

However, the opposite point of view to the possibility of existence of additional 
integrals is widespread in physical literature (see, for example, [I], where the ap- 
proach to justification of Gibbs distribution based on the formula (2.2) is called 
"speculative"). 

To show the possibilities of a new approach, let us consider a Hamiltonian system 
with n degrees of freedom. The Hamiltonian of this system has the form (1.4): 

H = H + eHx + o(e) , H = 'y]hj(x i ,y i ) , #i = . . . , x n , y u . . . , y n ) . 

(2.3) 

At e = the system splits into n independent subsystems with one degree of freedom. 
A system of connected pendula can be considered as an example. 

Let us assume that in each one- dimensional system with a Hamiltonian hi (1 < 
i < n) we can introduce the canonical action- angle variables Jj, ipi (mod27r) j^j. For 
example, for a system with the Hamiltonian h a from paragraph 1 such variables 
are defined in all phase space. In case of a pendulum the cylindrical phase space 
is divided by separatrices in three areas, in each area it is possible to introduce 
action-angle variables. In these variables 

h = h{Ii), l<i<n. (2.4) 

Let us assume that the functions (2.4) are continuous and monotonically increas- 
ing. The naturality of this supposition follows from the definition of an action vari- 
able as a normalized area on a phase plane contained inside a phase curve hi = const. 
Then the variable Jj will be a single- valued function of energy hi. 

The above mentioned assumptions are not fulfilled in the case, when the potential 
energy of a system has some local minimums. However, these conditions have a 
technical character, and probably they can be essentially weakened. 
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Thus, in action-angle variables the Hamiltonian (2.3) has a form 



H = J2 Hh) + eHtih, . . . , I n , (px, . . . , p n ) + o(e) . (2.5) 
i=i 

It is 27r-periodic with respect to each angular variable (pi, . . . , (p n . 

The density of probability distribution p for a system with the Hamiltonian (2.3) 
is a single- valued positive function of x, y, depending on parameter e. Let us assume 
that p is a function from class C 2 (it possesses the second continuous derivative on 
a set of In + 1 variables x, y, e). At small values of e this function can be written 

' ,S p = Po(x,y) +epi(x,y) + o(e) , (2.6) 

where po and p\ are functions from classes C 2 and C 1 respectively. 

Now let e tend to zero. Then po will be a density of probability distribution for 
a system with the Hamiltonian 

H = hi(x!, yt) + ... + h n (x n , y n ) . 

This system with n degrees of freedom splits into n independent systems with 
one degree of freedom. The main independence property denotes that the motion 
of each of these subsystems is uniquely determinated by any of its initial states. 

If we stick to the idea of statistical description of dynamical systems, it is nec- 
essary to introduce stationary densities of probability distribution 

Pi(an,2/i), . . .,p n (x n ,y n ) 

for each of one-dimensional systems. Taking into consideration the independence 
properties and using the probability product theorem, we obtain 

p =px...p n . (2.7) 

This equality is often called Gibbs hypothesis on thermodynamic equilibrium j2]. 

Remark 1. One should not think that any separation of variables results in 
independent systems with one degree of freedom. Let us show this fact by a simple 
example : 

\yl + xl) 



2H = y z 2 + x 



2 i Jl 
2 
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The variables x-i, y2 perform simple harmonic oscillations, the frequency of which 
is equal to the energy of oscillations of the first subsystem (described by canonical 
coordinates x\, y±). 

In the following, the Poincare set P [Zj plays an essential role. Let us decompose 
a disturbing function Hi in a multiple Fourier series: 

Hi = ^if (m) (/i, . . . , I n ) exp[i(mi(pi + ... + m n p n )} , m = (mi, . . . ,m n ) E Z n . 
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Let u)i = ^ (1 < i < n) be the frequencies of a nonperturbed problem; and 
assume uj = (ui, . . . , uj n ). By definition the set P consists of such points / = 
(Ii, . . . ,I n ) G M ra , for which there will be n — 1 linearly independent integer vec- 
tors a, a', . . . G Z n , such that 

1) (uj, a) = (uj, a!) — . . . — 0; 

2) ^ 0, F(«')(7) ^0,.... 

In a typical situation the set P fills in the range of values I G M n [Jj everywhere 
densely. 

Our basic result is the following: 

Theorem 1. Let us assume, that systems with the Hamiltonians hi, ... , h n are 
non- degenerate ^ ; the Poincare set is everywhere dense and the condition (2.7) 
is fulfilled. Then 

p = ce- pH0 {l + O{e)), (2.8) 

where c> 0, (3 are some constants. 

The constant c is inessential: the result of averaging with respect to the mea- 

SUre (2 - 8) (/} = jf Pdxdy 

J pdxdy 

does not depend on this constant. If s = the formula (2.8) gives the Gibbs 
canonical distribution (1.1). 

Remark 2. We stress the fact, that (2.8) is valid for the fixed number of degrees 
of freedom n > 2. A. A. Vlasov ' L 8l developed an approach to the derivation of Gibbs 
distribution, which would not use the analysis of interaction of subsystems at all and 
which is formally suitable for the case n = 1. This approach is based on the principle 
of "maximum statistical independence" , which seems to be an artificial supposition. 



Derivation of Gibbs distribution 



The proof of Theorem 1 is based on application of the Poincare method in the 
form indicated in [7j. 

Setting e = according to (2.6) we obtain that p is the first integral of a 
complete integrable Hamiltonian system with a Hamiltonian 

n 

k=i 

Since nonperturbed system is non-degenerate, i. e. 



det 



d 2 H 



dP 



n 

fe=i 



a*, 
dp ^ u 



then the function p expressed in action-angle variables /, <pmod 2tt depends only 
on 0. 
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Furthermore, since p is the first integral of the canonical system of differential 

equations with a Hamiltonian H, their Poisson bracket is equal to zero: {p, H} = 0. 

Let us differentiate this equality with respect to e and, then, let us assume e = 0. 

Since p and H are supposed to be the functions of a class C 2 the differentiations 

with respect to phase variables and to the parameter e are commutative. As a result 

we obtain the equality r „ 

{p ,H 1 } + {p 1 ,H } = 0, 

from which, with the help of the Fourier method and by a known method 7\, we 
can deduce that the functions p (/ 1 , . . . , J n ) and H = ^ hk(Ik) are dependent in 
all points of the Poincare set P. According to the supposition, this set is everywhere 
dense in a range of action variables /. Hence, the functions po and H are everywhere 
dependent by virtue of continuity. 

The variables Jj may be presented as single- valued functions of hi in accordance 
with the assumption made in Section 2 about properties of action-angle variables. 
Then p = Po(h\, . . . , h n ) and H = ^/ij. Since these functions are dependent, p 
is a smooth function of H . 

Indeed, the condition of dependence of po and H gives relations 

dp dp 



dhi dhj 



(3.1) 



which are valid for all values i, j. If we substitute h n = Hq — hi — . . . — h n -i in the 
expression for po, we shall obtain 

(n-l 
hi, ... , h n ~i, Hq — 2J hk 
fc=i 

However, this function does not actually depend on hi, ... , h n -i, since 

df dp dp 



dhi dhi dh 



, i < n 



n 



according to (3.1). 

The probability densities pi, . . . ,p n are the integrals of one-dimensional systems 
with Hamiltonians hi, . . . , h n . Hence, in the action- angle variables, Pi are smooth 
functions only of /j. But in this case pi is a single- valued differentiate function of 
the energy hi. 

Thus, the equality (2.7) can be presented in the following form : 

p (hi + . . . + h n ) = pi(hi) . ..p n {h n ) ■ 

Consecutively differentiating this relation with respect to hi, . . . ,h n and using the 
positiveness of functions pi, ■ ■ ■ ,p n , we obtain equalities 

^ = ... = ^ = -/3, 

Pi Pn 
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where (3 is some constant. Hence, 

Pi = c i e~ (3hl ' , Ci = const > , (3.2) 

and consequently 

po = ce~ pH ° , c = d . . . c n > . 

The theorem is proved. 

Remark 3. As we see from (3.2), the constant j3 is one and the same for all 
subsystems. It means that the weakly interacting subsystems are in the thermody- 
namic equilibrium ( and in the limit as e — > 0), since their temperatures T = are 
identical. Thus, when e — > 0, the statistical independence of subsystems is equiv- 
alent to their thermodynamic equilibrium. This is the physical sense of the Gibbs 
hypothesis. 

4 Analytical case 

In the applications, the Hamiltonian (3.2) is an analytical function of phase variables 
and the parameter e. In this case it is also natural to consider the probability 
density (2.6) as an analytical function with respect to x, y, e. 

The condition of everywhere- density of the Poincare set can be weakened: it is 
enough to require, that P would be a key set for a class of analytical functions. It 
means the following: if the analytical function f(I) is equal to zero in the points 
of P, then / = 0. The examples of nondense key sets may be found in In partic- 
ular, if the analytical functions are dependent in the points P, they are everywhere 
dependent. 

Theorem 2. Let us assume, that the systems with Hamiltonians h\, . . . ,h n are 
non- degenerate, the Poincare set is a key set for the class of analytical functions and 
the equality (2.7) is performed. Then the analytical density of probability distribution 
for the system with the Hamiltonian (2.3) has the form 

p = ce-P H [l + eg(H,e)}, (4.1) 

where c > 0, f3 ^ are some constants, g is an analytical function of the energy H 
and the parameter e. 

The function g can be represented as a power series in terms of e with coefficients 
depending on the energy H only. The summand eg in (4.1) is completely analogous 
to the Gram- harlier series in the theory of distribution of random variables which 
are unessentially different from the normal distributed ones (see, for example, [TO]). 

Proof. Let us prove Theorem 2. Expand the probability density in the power 
series of e: 

p = po + £pi + £ 2 P2 + ■■■ ■ 
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On account of non-degeneracy of the unperturbed system, p is an analytical function 
of action variables / only. Since po an d H are dependent at points / e P, and P is the 
key set, po and Hq are dependent everywhere. Therefore, taking into consideration 
the Gibbs hypothesis (2.7), 

po = ce" 13110 ; c,(3 — const, 

(see Section 3). Since p > 0, we have c > 0. The constant (3 is not equal to zero, 
otherwise po is not the density of probabilistic measure (the volume of the whole 
phase space is infinite). 

Further, the analytical function 



ce 



-f3H 



is an integral of Hamilton equations with the Hamiltonian (2.7). Let us expand this 
function in power series of e: 1 + eGo + e 2 Gi + . . . . It is obvious that the series 

G + ed + . . . (4.2) 

is an integral of the same system. With the help of the method used in Section 3 
one can prove that Go is an analytical function of H : Go = go(H ). It is obvious 
that the power series 

[G . +eG , + ..- a ,( g )] = 

e 

is again a first integral. Thus, F Q is an analytical function of H : F = gi(H ). 
Extending infinitely this procedure, we obtain that the series (4.2) has actually the 
form: 

g (H) + £ 9l (H) + . . . . 

Denoting this function by g(H,e), we obtain the desired formula (4.1). ■ 
Remark 4. Let us introduce K(H,e) by setting 

p = ce"^, (4.3) 

where the density p is given by (4.1). The function K can be expanded in power 
series K + eK 1 + . . ., where 

K = H , K 1 = Hl - 9 -^,.... 

Hamilton equations with the Hamiltonians H and K possess the same trajectories, 
however, times of motion along these trajectories are different. The formula (4.3) 
implies the Gibbs canonical distribution for a slightly changed Hamiltonian system. 
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5 Application to the system of weakly connected 
pendula 

Let us consider n identical mathematical pendula of mass m and length I, con- 
secutively connected with each other by elastic springs with small rigidity x. For 
the simplicity we assume that the fixed points of pendula coincide. This system 
with n degrees of freedom is described by canonical differential equations with the 
Hamiltonian Hq + eH\ , where 

n-l 

^COs(x; - X i+1 ) , 
i=l 

e = — is a small parameter. The transition to angle-action variables is carried 
out for the pendulum with the help of elliptic functions (see, for example, It 
is possible to show that the Fourier series of perturbation function with respect to 
angle variables (p\, . . . ,(p n mod 2n has the form: 

tt _ \ r L / j \ 2i(m 1 <p 1 +m 2 tp2) i i \ r /, I T \ e ,2i(m n - 1 (p n -i+m n <Pn) 

1A 1 — / j n m 1 ,m2\ 1 l) c i ■ ■ ■ i / j ll -m n ^ 1 ,m n \ 1 n)t 

rai,m2 m n -i,m n 

The summation is taken over all integer mi, . . . ,m n from — oo to +00. The coeffi- 
cients in this expansion can be expressed explicitly with the help of known expansions 
of elliptic functions sn 2 , cn 2 and sncn in Fourier series [TI]. All of them are not 
equal to zero. 

The Poincare set P is defined in this problem as the set of points I = . . . , I n ), 
satisfying n — 1 equations 

mia;i(ii) + m 2 u; 2 (7 2 ) = . . . = m n _i^ n _ 1 (J n _i) + m n u n {I n ) = , 

where either m!m 2 . . . m n _i ^ or m 2 m 3 . . . m n 7^ 0. If one of the two lat- 
ter conditions is fulfilled, the vectors (mi, m 2 , 0, . . . , 0), (0, m 2 , m.3, 0, . . . , 0), 
(0, 0, . . . , m n _i, m n ) are linearly independent. It is possible to show that the set P, 
which consists of infinitely many curves, fills in the domain of definition of action 
variables {Ji > 0, . . . , I n > 0} everywhere densely. 

Thus, the statements of Theorems 1 and 2 are valid for the chain of connected 
pendula. If one of the springs is taken away, the system splits into two disconnected 
chains, the Hamiltonians of which are the first integrals of the total system. In this 
case the density of probability distribution is not the function of the total energy of 
the system only and thus is not subjected to the Gibbs distribution. 



2m/ 2 



mgl cos Xi 
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6 Thermodynamics of mechanical systems 



Let us consider a Hamiltonian system with the Hamiltonian H, where the den- 
sity of probability distribution is defined by the Gibbs formula (1.1). According 
to theorem 1, it can be a system with n degrees of freedom, composed of inde- 
pendent one-dimensional subsystems, or one of these subsystems. We shall show, 
after Gibbs |T], that this mechanical system can be naturally connected with some 
thermodynamic system. 

Let us assume that the Hamiltonian H depends not only on canonical variables 
x, y, but also on several parameters Ai,...,A m . One can take, for example, the 
length of the pendulum as such a parameter. The parameters A can be considered 
as generalized coordinates of some system with n + m degrees of freedom and the 
constancy of A as application of holonomic constraints. Thus, let H*(x, y, A, /i) be 
the Hamiltonian of a system with n + m degrees of freedom, /zi, . . . , /i m be canonical 
momenta conjugated with additional coordinates Ai, . . . , A m . The dynamics of the 
extended system is described by the canonical equations 

oy ox o[i OA 

Let us impose m independent relations Ai, . . . , A m = const on this system. There- 
fore, Aj = 0, and from the equations 

oti* oH* , . 

—JL = ... = — ^ = (6.2) 

dill OH m 

one can derive momenta fi as functions of x, y and constant parameters A. The 
sufficient condition of solvability of the equations (6.2) with respect to \i is reduced 
to the inequality 

d 2 H* 



det 



^0. 



d/iidfij 

It is fulfilled automatically if is a positively defined quadratic form of n + m 
momenta y, fi. 

Substituting the obtained expressions for \i in the last set of equations (6.1), 
we obtain additional relations on canonical variables x, y, which, of course, are not 
fulfilled in general case. Therefore, it is necessary to introduce additional forces, the 
constraint reactions Ri, . . . , R m , and replace the latter equation of (6.1) by 

H=-^ + R. (6.3) 
Let \i = fi(x,y, A) be a solution of the algebraic system (6.2). We assume that 

H(x, y, A) = H*(x, y, A, fj,(x, y, A)). 
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On account of (6.2), 

dH _ dH* dH _ dH* dH _ dH* 
dx dx ' dy dy ' dX dX 

Consequently, for constant values of A the variables x, y change according to Hamil- 
ton equations with the Hamiltonian H, and the equation (6.3) can be replaced by 

p = -°— + R. (6.4) 

Let x(t,x ,yo), y(t,x ,y ) be solutions of canonical equations with the Hamilto- 
nian H. We assume that for every such solution fi(t)/t — > as t — > +oo. Then the 
time-average of (i is equal to zero. This assumption is automatically fulfilled if the 
configuration space {x} is compact: on account of existence of energy integral the 
function p{t) is bounded. 

At first let us average both parts of the equality (6.4) with respect to time 
and then with respect to the measure p(x ,y )dx dy , where p is given by (1.1). 
According to the Birkhoff-Khinchin ergodic theorem ^2] the obtained relation is 
equivalent to the following one: 

<*> = (§). (6,) 

Here ( ) is the mean with respect to the Gibbs measure (1.1). 

Remark 5. Usually J3|,^ the relation (6.5) is derived from the simplified relation 
0/(6.4), which lacks the derivative fi. 

The relation (1.2) defines the "internal" energy E as a function of parameters 
Ai, . . . , A m and f3 (recall that f3~ l = kT). Let us set A = —(R), i. e., the phase mean 
values of the above introduced constraint reactions (with the opposite sign). The 
relations 

Ai = fi(Xi, . . . , X m ,(3) , l<i<m (6.6) 

are usually called in thermodynamics the equations of state. Assignment of func- 
tions E and Aj is included in the definition of thermodynamic system. However, 
these functions can not be arbitrary, since the first and the second principles of 
thermodynamics must be fulfilled. 

Let us introduce the statistical integral 



Z{\,/3) = J e~ pH dxdy. 



One can verify the validity of the equalities [3 E] 



IdlnZ . . „ d\nZ , m . 
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In thermodynamics the main role plays the differential 1-form of heat gain 

uj = dE + AjdXj . 

Motivations for this definition can be found in [311] ■ According to the first equality of 
(6.7) the form to is the exact differential for fixed values of f3 (or absolute temperature 
T). This is the first principle of thermodynamics. Then, taking into account (6.7) 
we obtain 

f3uj = (3dE + J2 A i dX i = d^E) - Edf3 + J2 A i dX i = d (P E + 

Thus, the form of heat gain possesses the integrating factor (3 = l/(kT). This is the 
second principle of thermodynamics. The function S = f3E + In Z is the entropy of 
thermodynamic system. 

Finally, let us give an illustrative example. We consider a motion of a point of 
mass m, attached to the end of unstretchable thread of length I. Let x be the angle 
of rotation of the thread, y be the conjugate canonical momentum. If active forces 
are not applied to the point, its dynamics is described by canonical equations with 
the Hamiltonian 

2ml 2 ' 

Let us take the length of the thread I as the parameter A. The statistical integral is 
equal to 

Z = (2nf 2 m l /H/3- 1 / 2 . 

Let p be a force, corresponding to the parameter I (the thread tension). Then, 
according to (6.7) the equation of state (6.6) has the form p = l/(/3l) or pi = kT. 
It is similar to the ideal gas law. The second equation (6.7) gives the relation for 
the internal energy E = 1/(2/3) or E = kT/2. From these relations we obtain the 
formula p = 2E/1, which, however, was valid before the application of averaging 
procedure. 

The work is partly supported by RFBR (No. 99-01-01096) and INTAS (No. 96- 
0793). 
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